Forcing anomalous scaling on demographic fluctuations 
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We discuss the conditions under which a population of anomalously diffusing individuals can be 
characterized by demographic fluctuations that are anomalously scaling themselves. Two examples 
are provided in the case of individuals migrating by Gaussian diffusion, and by a sequence of Levy 
flights. 
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I. INTRODUCTION 

Simple birth-death models, with individual migration 
described by a diffusion process, provide a simple illustra- 
tion of how demographic stochasticity at the microscopic 
level translates into the phenomenon of spatial clustering 
[1, 2]. The original setting of the model involved Brown- 
ian walkers ( "Brownian bugs" , in the terminology of [3] ) , 
but the generalization to the case of anomalously diffus- 
ing individuals was considered as well [4, 5] . Contrary to 
what could have been expected, clustering did not appear 
to be sensitive to the character of the diffusion process. 
What appeared to be effective in modifying the popula- 
tion dynamics, instead, was the possibility that memory 
of the trajectories be transferred from individuals to their 
descendants [5]. 

At least in the case of a migration dynamics of the 
continuous time random walk (CTRW) type [6, 7], mem- 
ory transfer between generations was insufficient to pro- 
duce any scaling dependence on the Hurst exponent of 
the diffusion in the population dynamics [5]. What was 
observed, was collapse on a Galton- Watson like behav- 
ior [8] . The question remained open on whether different 
models would lead to the same result. 

What we intend to analyze in the present report is pre- 
cisely under what conditions a dependence between the 
scaling of the local population fluctuations, and that of 
the diffusion process, can exist. We shall consider two 
cases: that of anomalous diffusion produced by a Gaus- 
sian process, and that of diffusion produced by a sequence 
of Levy flights [9, 10]. These are situation of impor- 
tance both conceptually and from the point of view of 
applications. Physically relevant examples of anomalous 
Gaussian diffusion include particles migrating with a ve- 
locity that is solution of a generalized Langevin equation 
(GLE) [11], and, of course, the fractional Brownian mo- 
tion (FBM) [12]. Levy flights constitute an important 
ingredient in modelling the spreading of epidemies [13], 
biologic searching strategies [14], as well as bacteria dy- 
namics [15]. 

In the case of Gaussian diffusion, we shall see that 
"anomalous" scaling of the population fluctuation can 
indeed be produced, and it is the result of a delicate 
interplay between Gaussian statistics, memory transfer 
between generations, and the non-Markovian nature of 



the process. In the case of Levy flights, the Markovian 
nature of the diffusion process makes any memory trans- 
fer between generations irrelevant to the dynamics, and 
the scaling of the local population fluctuations will ap- 
pear to be determined automatically by the stable law of 
the Levy process. 



II. GAUSSIAN DIFFUSION 

Let us consider the simple situation of a population 
of non-interacting individuals, with identical birth and 
death rates Ts = Td = T (one offspring for each birth 
event) to insure stationarity. We assume births and 
deaths to be Markovian and focus on the effect that 
anomalous diffusion of the individuals may have at the 
population level. Let us limit the analysis, for the mo- 
ment, to a one-dimensional situation. 

We consider first the case of a Gaussian, anomalous dif- 
fusion process. Contrary to the case of normal diffusion, 
in which the condition of linear scaling of the mean square 
displacement (assuming Gaussian statistics) determines 
uniquely the process, full knowledge of the displacement 
correlation is necessary in the anomalous case. 

Let us consider a population that initially is uniformly 
distributed in space, with pi(x, 0) = uq. Equality of 
the birth and death rates guarantees that the population 
remains on the average uniformly distributed, with the 
same population density p\{x,t) = (h{x,t)) — no [we 
have indicated with n(x, t) the instantaneous value of 
the population density, including fluctuations] . To study 
fluctuations, we must consider the higher moments of 
h(x,t), in particular the two-bug correlation 



P2(xi,x 2 ;t) = (n(xi,t)n(x 2 ,t)), 



(1) 



or, more precisely, its connected part p2c(xi,X2',t) — 
P2{x\,x 2 \t) - n§. 

A representation of the population dynamics is ob- 
tained considering the family trees of its members, as 
illustrated in Fig. 1. A quantity such as the cor- 
relation p2{x\,X2]t) is obtained therefore by summing 
over the tree branches passing through [x\,t) and (x2, t) 
and averaging over the family tree configurations. The 
family trees responsible for the connected component 
P2c(xi,X2]t) are those represented in case a of Fig. 1, 
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FIG. 1: Family tree contributing to the connected (a) and 
disconnected (6) part of the two-bug correlation p2(xi, xi; t). 
[see Eq. (1)]. The correlation is obtained counting pairs of 
branches passing through the sampling points at a?i,2- 



in which the individuals at (2:1,2, t) share a common an- 
cestor. Clearly, cutting the family tree of case a above 
the branching at (z,t) and neglecting all the previous 
history, would lead to a disconnected tree identical to 
the one represented in b. It is easy to be convinced that 
the same result would be obtained by reinitializing from 
scratch the migration strategies of the individuals at the 
branching. A complete reinitialization of migration at 
each birth event, would lead in fact to Markovianization 
of the population dynamics at times longer than the in- 
dividual lifetime with the displacements of different 
individuals over a lifetime becoming independent. (In the 
case of migration by a CTRW, this lead to a clustering 
dynamics in the population identical to that of Brownian 
bugs [16]). 

In order to obtain a non-trivial dynamics, the new- 
born individuals must preserve some memory of the tra- 
jectories followed by their parents. We shall consider the 
situation in which individuals along the same family line 
have full memory of the trajectory followed by their an- 
cestors. 

Let us indicate by G{x — z,t — r) = pi(x, t\z, r) the 
mean population density at (x,t) given presence of an 
ancestor at (z, r), with r < t, and with pi(x, t\z, r; w, 0), 
the same quantity conditioned to presence of a common 
ancestor at (w, 0). From equality of the birth and death 
rates, the total number of individuals accounted for by 
pi(x,t\.) is one and these functions are normalized to 
one and coincide with the conditional PDF's (probability 
density functions) for the position a single walker in the 
absence of birth and death effects. We can write: 



0) 



p2c{xi,x 2 ;t) = 2T J dr J dz J dw p 1 (x 1 ,t\z,r;w, 

x pi(x2,t\z,r; w, 0)G(z — w, r) 

x piKO). (2) 

We can read Eq. (2) from the family tree a in Fig. 1: 

• The factor G(z — w,t)p\(w,0) is the individual 
density at (z,t), originating from an ancestor 
at (w,0), that is available to generate offspring; 



TG(z — w, r)pi(w, 0) is the corresponding total rate 
of birth. 

• The factors p\ (#1,2, t\z, r; w, 0) enforce memory of 
the position of the common ancestor at (w, 0) on 
the offspring at (z,t), until it (or one of its descen- 
dants) reaches 2:1,2. 

• The factor 2 on the right side accounts for the num- 
ber of ways in which parent and offspring at (z, r) 
(or their descendants) may distribute between x\^. 

Gaussian statistics makes determination of the condi- 
tional PDF's pi(x,t\.) straightforward [19]: 



G(z — w, t) 



1 
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■ exp 
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pi(x,t\z,T;w,0) 



x exp 



1 



" (2n) 1 / 2 a(t\T;z-w) 
\x — p.{t\r, z — w)\ 2 
2a 2 (t\T;z -w) 



(4) 



where the conditional mean and mean square displace- 
ments n(t\.) and cr 2 (i|.) are given by 



/i(t|r, z - w) = w + ^ v " (z - w) 



and 



2/,, x _ 2/, \ 2/.n (yitMr)) 

a (i r; z - w) = a (i; t) = a (t) ^— — . 



(5) 



(6) 



with y(t) = x(t) - a;(0) and <r 2 (t) = (y 2 (t)). 

Substituting Eqs. (3-6) into Eq. (2), and using 
Pi(x, 0) = no, we obtain the result 



P2c(xi,x 2 ;t) 



Trip 
2^ 



r dr 
7T Jo °(^ T 



exp 



(xi - x 2 ) 2 



Aa 2 {t,r) 



•(7) 



For i » t we see that a(t;r) — > cr(t). The Markovian 
case is obtained substituting a(t; r) — > a(t — r), for t, r 
generic, in Eq. (7), that becomes in this way the solution 
of a heat equation with constant Dirac delta forcing at 
x\ = x 2 [16]. 

The scaling of the square amplitude fluctuation 
p 2c (x,x;t) can be calculated explicitly if the correlation 
profile entering Eq. (6) is known. In the case of an FBM 
[12]: 



t 



(y(t)y(r))=(K H /2)(t 2H +t 
and we get 

p 2c (x,x;t) = CTnoK^t 1 ^ 11 

with 

1 f 1 u H du 



l2H\ 



(8) 
(9) 



C = 



2V^Jo ^u 2H - [1 + u 2H - (1 - u) 2H ] 2 /2 



III. LEVY FLIGHTS 
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FIG. 2: Fluctuation build-up in the case the individuals mi- 
grate by Gaussian diffusion (heavy line) and by Levy flights 
(thin line). The slope t 1 '" , for H = 1//3 = 0.75, is shown 
for reference. In both cases N = 10 5 individuals in a peri- 
odic domain were considered. In the time units considered: 
the fastest mode in the FBM has relaxation time one; the 
discretization of the Levy flights is At = 0.25; Y = 0.1. 



The general condition <7(t;r) — > <r(t) for t/r large, leads 
us to expect that the same scaling of Eq. (8) be recov- 
ered also for other Gaussian, anomalous diffusion pro- 
cesses. An example is illustrated in Fig. 2, in the case 
of a supcrdiffusive process generated with the algorithm 
described in [17], in which the velocity of the individu- 
als is obtained as a superposition of Ornstein-Uhlenbeck 
processes. 

From Eq. (7) we could also obtain the scaling of 
the correlation length of the fluctuations, A c : X 2 = 
[/ dy p 2c (x, y; t)}- 1 J dy (x - y) 2 p 2c (x, y; t). Wc obtain: 

Ac(t) ~ 4 /2 i ff . (10) 

We thus see that the slow-down of the fluctuation build- 
up in Eq. (9), with respect to a case without migra- 
tion: p2 C {x,x;t) = 2n Tt (Galton- Watson dynamics [8]), 
is precisely the smearing of the fluctuations that is pro- 
duced by anomalous diffusion. We notice that clustering 
is possible (in one dimension) only as long as H < 1; thus 
superdiffusive processes that are faster than ballistic are 
excluded. 



In this case we do not have to bother with memory 
transfer among generations, as the process is Markovian. 
The process can be obtained as the continuous limit of 
a dynamics in which the individuals, at time intervals 
At, carry on independent jumps distributed with a PDF 
characterized by power-law tails: 

G(y,At)~ M- 1 ^, 0</3<2 (11) 

(Levy flights). In Fourier space, this corresponds to the 
asymptotic behavior for small wavenumbers: 

G k {At)~l-c\kfAt, (12) 
where c ~ (Ax) 13 / At, and Aa; is a microscale below which 
the scaling in Eq. (11) ceases to hold. The slow decay in 
y of the PDF, Eq. (11), implies divergence of the second 
moment of the displacement (\y{t)\ 2 ) — oo even for a 
single jump, and the PDF G(y, t) can be shown to tend 
to a Levy distribution; in Fourier space [6]: 

G k (t) = exp{-c\kft). (13) 

Equation (13) tells us that, although G(y,t) does not 
have a second moment, we can identify a characteristic 
scale: 

y(t) = x(t)-x(0)~(ct)^, (14) 

with the exponent 1 / (3 playing a role analogous to that of 
the Hurst exponent for a diffusive process (more precisely, 
for a superdiffusive process, as 1/(3 > 1/2). 

As in the previous section, let us assume that the indi- 
viduals undergo processes of birth and death with equal 
rate T, and that the individuals are distributed at time 
t = with uniform density no. As in the previous case, 
it is easy to see that the mean density remains constant 
pi(x,t) = no, but fluctuations of increasing amplitude, 
accounted for by the connected correlation p 2c (x\,x 2 ;t), 
are generated. 

The Markovian nature of the process allows us to de- 
rive a version of the evolution equation for p 2c , Eq. (2), 
that is local in time: 



P2c(xi,x 2 ;t + At) = J dyi J dy 2 G(x± - yi, At)G(x 2 - y 2 , At)p 2c (y\,y 2 ;t) + 2Tn At5(xi - x 2 ). (15) 

I 

Fourier transforming in xi j2 and defining p 2c .k 1 k 2 {t) = can t nus use Eq. (12), and Eq. (16) becomes, taking the 
27r<5(fci + k 2 )Ck 1 (t), we obtain continuous limit: 

C k (t + At) = G 2 2k (At)C k (t) + 2Tn a At. (16) C k + 2c\kfC k = 2Tn a , (17) 



For t 3> At, the space scale contributing in Eq. (16) are that can be seen as a form of forced fractional heat equa- 
those corresponding to the asymptotics in Eq. (11). We tion [18]. Solution by Laplace transform in time gives 
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us: 

Ckz = Z { Z 2 +2c\kn (18) 

Inverse Fourier transform at zero space separation gives 
then the result 

CM = J ^c kz = BTn c-ve z ve- 2 , 

where B = j* + °°(l + 2h)- 1 h 1 ? fj - 1 dh, and this corre- 
sponds to the long-time asymptotics 

p 2c (x,x;t)^Tn a c- 1 ^t 1 - 1 ^. (19) 

A realization of such a process is shown in Fig. 2. Again, 
fast processes for which (3 < 1 do not lead to clustering. 

In spite of the fact that the dynamics considered in 
this and the previous section are completely different, as 
clear from Eqs. (9) and (19), the fluctuation build-up 
occurs in the two cases in the same way. Looking at Eq. 
(18) and comparing with Eqs. (7) and (10), we find in 
fact the same mechanism of smearing of fluctuations at 
scale x(t) cx t 1 ^ [Eq. (14), Levy flights] or X c (t) cx t H 
[Eq. (10), Gaussian diffusion]. 

IV. CONCLUSION 

We have provided two separate examples of how the 
scaling in the build-up of demographic fluctuations in a 
birth-death model, can be made dependent on the kind 
of diffusion utilized for migration. We have seen that 
the basic ingredient for such dependence is the smear- 
ing of the fluctuations by the diffusion process, provided 
the inclusion of demography does not affect migration at 
population scale. With this we intend the fact that the 
dispersion of a group of individuals has the same scaling 
in time independently of whether demography is taken 
into account or not. 

As in the case of Brownian bugs, working in higher 
dimension D > 1, would decrease the build-up exponent 



in proportion; namely, in place of Eq. (9) or (19), it is 
possible to see that we would have: p2 C {x,x;t) cx t x ~ DH 
or p2c(x, x; t) cx t 1 ^ ^ in the two cases. We thus see that 
superdiffusion would lead to clustering only for D = 1. 
However, in contrast with the case of Brownian bugs, 
clustering may arise in D > 1 in the case of sufficiently 
slow subdiffusion (at least if diffusion is produced by a 
Gaussian process). 

Summarizing, and including the results in [5] as well, 
the various dynamical mechanism leading to anomalous 
diffusion, and the different ways in which memory is 
transferred between different generations of individuals, 
produce the following results: 

• Gaussian anomalous diffusion (FBM, individuals 
moving with velocity that is solution of a GLE, and 
others): if offsprings share memory of the trajecto- 
ries with their parents, demographic fluctuations 
will scale anomalously, as described by Eq. (9). 
Otherwise, the case of Brownian bugs is recovered. 

• Levy flights: demographic fluctuations scale 
anomalously, as described by Eq. (19). 

• CTRW: in the absence of memory transfer between 
generations, the case of Brownian bugs is recovered. 
If the offsprings share their escape time (the only 
thing that they can share) with their parents, the 
dynamics falls back on that of a Galton- Watson 
process. 

• A spatial assembly of random traps: the same be- 
havior of the CTRW with memory transfer between 
generation is obtained; the dynamics falls back on 
that of a Galton- Watson process. 
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